#include <complex>
inline std::pair< std::complex<double>, std::complex<double> > Quad(double a, double b, double c)
{
    std::complex <double> square = (b*b * -4*(a*c));
    return {(-b + std::sqrt(square)/(2*a)), (-b - std::sqrt(square)/(2*a))};
}
